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liminary CDF upper limit on double diffractive exclusive x° c production 
cross section. 
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1 Introduction 

The study of the double pomeron exchange (DPE) production processes is interest- 
ing in its own right. It is an ideal way to improve our understanding of diffractive 
processes and the dynamics of the pomeron exchange. 

However, the great interest in such reactions is caused by the possibility of the 
DPE processes to be one of the main mechanisms leading to Higgs boson production 
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[1-10] within a very clean experimental environment. 

In the present Letter we are particularly interested in the exclusive and central 
inclusive (central inelastic) DPE production of heavy quarkonium states Xc and Xb 
[11-13], see also [14, 15]. In the exclusive DPE event the central object x is produced 
alone, separated from the outgoing hadrons by rapidity gaps: 

pp — > p + gap + x + g a P + P- 

In the central inclusive DPE event an additional radiation accompanying the central 
object is allowed. 

The basis for our considerations is the Bialas-Landshoff model for central in- 
clusive double diffractive Higgs boson production [2]. We showed that the Bialas- 
Landshoff model and its exclusive extension [10] give satisfactory description of the 
DPE central inclusive and exclusive dijet cross sections [10, 16]. 

In this Letter we show that the exclusive extension of the Bialas-Landshoff model 
is consistent with the preliminary CDF upper limit on double diffractive exclusive 
X® production cross section [17]. 
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2 Central inclusive % meson production 



In the Bialas-Landshoff approach pomeron exchange corresponds to the exchange 
of a pair of non-perturbative gluons which takes place between a pair of colliding 
quarks [18]. Thus, the process of x meson production is described by a sum of 
the four (reggeized) diagrams shown in Fig. 1 where dashed lines correspond to 
the non-perturbative gluons. The x coupling is taken to be through a c-quark and 
6-quark loop for Xc and Xb respectively. 




Figure 1: Four diagrams contributing to the amplitude of the process of x meson 
production by double pomeron exchange. The dashed lines represent the exchange 
of the non-perturbative gluons. The x coupling is taken to be through a c-quark 
and 6-quark loop for Xc and Xb respectively. 

Our calculation follows closely that of Ref. [2] where the DPE contribution to 
the central inclusive Higgs boson production is calculated. 

First, one calculates the (non-reggeized) production amplitude in the forward 
direction i.e. for vanishing transverse momenta of the produced x meson and of the 
final hadrons. In this case the sum of the diagrams of Fig. 1 can be approximately 
replaced by the s-channel discontinuity of the first one [2], shown in Fig. 2. 

In the second step one introduces phenomenologically the effects of reggeization, 
transverse momentum dependence and the gap survival so that the full amplitude 
becomes: 1 

/ \ a(t 2 )-l / \ 

M=M I— J I -J F(t 1 )F(t 2 )exp(/3(t 1 + t 2 )) 1 S gap . (1) 

Here A4 is the amplitude in the forward scattering limit given by the diagram shown 
in Fig. 2. a (t) = l+e+a't is the pomeron Regge trajectory with e w 0.08, a' = 0.25 

1 For DPE central inclusive %c and \b mesons production we can neglect the additional gap 
spoiling effect i.e. the Sudakov effect. It is justified by the relatively small masses of the produced 
mesons. 
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Figure 2: Putting the inner quark lines "on-shell" is equivalent to calculating the 
amplitude for \ meson production by double pomeron exchange in the forward 
scattering limit i.e. t± — t 2 — 0. 

GeV~ 2 . s = (P1+P2) 2 , si = (h + P) 2 , s 2 = (k 2 + P) 2 , h = ( Pl -h) 2 , t 2 = (p 2 -k 2 ) 2 
where pi, p 2 , hi, k 2 and P are defined in Fig. 2. F (t) = exp(At) is the nucleon 
form-factor with A = 2 GeV~ 2 . The phenomenological factor exp (f3 (ti + t 2 )) with 
(3 = 1 GeV -2 takes into account the effect of the momentum transfer dependence 
of the non-perturbative gluon propagator given by (p 2 is the Lorentz square of the 
momentum carried by the non-perturbative gluon): 2 

D{p 2 ) =D exp{p 2 /r 2 ), (2) 

with t = 1 GeV and D G 2 r = 30 GeV -1 [2] where G is the scale of the process 
independent non-perturbative quark gluon coupling. 

The factor S gaiP takes the gap survival effect into account i.e. the probability 
(5gap) of the gaps not to be populated by secondaries produced in the soft rescat- 
tering. It is not a universal number but it depends on the initial energy and the 
particular final state. Theoretical predictions of the gap survival factor S 2 ap can be 
found in Ref. [19]. In our calculations, following [10,16], we take for the Tevatron 
(LHC) energy S 2 g J (G 2 /4tt) 2 = 0.6 (0.25). 

Before we present the details of our calculation one point must be emphasized. 
Here we calculate in a model dependent way the non-perturbative part of the whole 
amplitude. The contribution coming from small distances, being another part of 
the whole amplitude, is beyond our approach. However, for the case of x nieson 
production the main part of the cross section may has a non-perturbative origin. 

Following the calculation presented in Ref. [2] we find Aio for colliding hadrons 3 

2 As was stated in Ref. [2] there is no reason to believe that the true form of D is as simple as 
this. We hope it is not a serious objection to our model. 

3 The calculation of Mo for colliding hadrons is performed in two steps. First, we calculate the 
diagram presented in Fig. 2 for quark-quark scattering. Next, we multiply by a factor 3 2 to take 
the presence of three quarks in each (anti)proton into account. 
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to be in the form: 4 

Mo = J d 2 Q T ^^exp(-3g 2 T /r 2 ). (3) 

Here Q T is the transverse momentum carried by each of the three gluons. V^ v is 
the gg — > x J vertex depending on the polarization J of the \ J meson state. It was 
shown [13] that the DPE contribution to x 1 an d X 2 production in the forward scat- 
tering limit is vanishing (either perturbative or non-perturbative two gluon exchange 
models). 

It turns out that for J = we have the following simple result [2, 11, 12]: 

rfVft* = $-A, (4) 

x° 

where A is expressed by the mass M x oand the width T x oof the x° meson through 
the relation: 

^° = 2!4^" (5) 

Substituting (3) with (4) and (5) to the full amplitude (1) and then performing 
the appropriate calculations [16] we find the differential cross section da/dy for DPE 
central inclusive x° nieson production to be in the form: 

do _( s CR 2 
dy ~ \M[°) ((A + /?)/«' + ^ iV~s/M x »}f-yi' 
Here y, ln(M x o / '{Sf^y/s)) ^ y < ln(5™^v^/M x o), is a rapidity of the produced x° 
where 5i and 5 2 are defined as 5± j2 = 1 — ^1,2/^1,2 k 2 , Pi, p 2 , are defined in Fig. 
2). In the following we take c^ ax = Sf** = 5 = 0.1. 

The factors C and R are defined as: 

° = (12^Mj^ ( DoG T) T (GW' ( } 

R = 9 J dQ 2 Q 2 exp(-3Q 2 T ) = 1. (8) 

In the expression (6) R 2 reflects the structure of the loop integral and is shown 
explicitly for the reason which will become clear in the following. 

Performing integration over the whole range of rapidity we obtain the following 
result for the central inclusive total cross section: 



2c 



CR 2 , (X + P)/a' + \n(Ss/M 2 



M 2 J (\ + fi)/a' + In [sfs/M x o\ ^ y (A + /3)/a' - In 5 ' (9) 

This completes the calculation of the cross section. 

4 This formula is only valid in the limit of 61,2 << 1 where #1,2 = 1 — ki^/pi, 2, and for small 
momentum transfer between initial and final quarks. 
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3 Sudakov Factor 



The calculation presented in the previous section, based on the original Bialas- 
Landshoff model, is a central inclusive one, i.e. the QCD radiation accompanying 
the produced object is allowed. 5 Thus, to describe the exclusive processes where the 
central object is produced alone one has to forbid this radiation. To this end we 
shall include the Sudakov survival factor T(Q T ,fi) [5] inside the loop integral over 
Q T (8). The Sudakov factor T(Q T ,fi) is the survival probability that a gluon with 
transverse momentum Q T remains untouched in the evolution up to the hard scale 
(i = M/2. 

Naturally a question of internal consistency arises. Namely, the Sudakov factor 
uses perturbative gluons whilst in our calculations of the Born amplitude (3) we 
used non-perturbative gluons. We hope however that taking the Sudakov factor 
in the loop integral into account we obtain an approximate insight into exclusive 
processes. Moreover, it is shown [10] that such approach leads to the satisfactory 
description of the DPE exclusive dijet production cross sections [20,21]. 

We take the function T(Q T ,[i) to be in the form [5]: 



T(g ^H-4 * if/ 



dz . (10) 



Here A = |fc T | / (// + |fc T |), Pgg{ z ) an d P qg {z) (we take q = u,d,s,u,d,s) are the 
GLAP spitting functions. a s is the strong coupling constant. 6 Taking into account 
the leading-order contributions [22] to the GLAP splitting functions: 

P ga (z) = 6 [z/(l -z) + (l- z)/z + z (1 - z) + S (1 - z) (11/2 - n/ /3)] , 
P q9 (z)=[z 2 + (l-z) 2 }/2, (11) 

we obtain: 

-l-A 



l-A 



/ zPgg (z) dz = -11/2 + 12A - 9 A 2 + 4A 3 - 3A 4 /2 - 6 In A, 
Jo 

P qg (z) dz = 1/3 - A/2 + A 2 /2 - A 3 /3. (12) 
Now to describe the exclusive processes we use the formulas (6,9) with R 2 = 1 



I 



5 Let us notice that the central inclusive process contains exclusive process. It implies that the 
relation cr c .inci. > cr CX ci. holds. 

6 In the following we take a s at one loop accuracy i.e. a s (q 2 ) — (47r//3 ) (l/ In (g 2 /A 2 )) with 
[3q = 9 and A = 200 MeV. Below q = 0.8 GeV we freeze a s to be 0.5. 
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replaced by R 2 (//) where R (//) is defined as: 7 

R (/i) = 9 jT dQ 2 Q 2 exp (-3g 2 T ) T (Q T , /i) . (13) 

4 Exclusive x meson production. The CDF result 

Now we are ready to give our predictions for DPE exclusive Xc an d Xb mesons 
production at the Tevatron and the LHC energies. The masses and widths of the 
X° c and xl are taken as follows: M x o = 3.4 GeV, M x o = 9.8 GeV, T x o = 15 MeV, 
r x o = 2.15 MeV. Taking into account (13) we find the effective Sudakov suppression 
of the cross section to be (/i = M/2): 

R 2 ( X ° C ) = 0.3, 

R\xl) = 0.07. (14) 

In Table 1 the predictions for the Tevatron and the LHC energies are shown. As 
was discussed earlier we take 5 = 0.1 and assume S 2 ap / (G 2 /Att) 2 to be 0.6 (0.25) for 
the Tevatron (LHC) energy. 





da/dy(y = 0) 


a 


2 TeV 
Xc 14 TeV 


45 nb 
30 nb 


370 nb 
350 nb 


2 TeV 
Xb 14 TeV 


0.05 nb 
0.03 nb 


0.3 nb 
0.3 nb 



Table 1: Our results for DPE exclusive x% Xb niesons production cross section for 
the Tevatron and the LHC energies. 

As can be seen from Table 1 the obtained results for the LHC energy are com- 
parable or even smaller than those for the Tevatron energy. It is due to the faster 
decrease of the rapidity gap survival factor S 2 ap with increasing energy than the in- 
crease of our "soft" energy dependence s 2t . It should be noted, however, that the s 
dependence is the most serious uncertainty of our approach and the results presented 
in Table 1 should be regarded only as an order of magnitude estimates. Bearing in 
mind this uncertainty our results seem to be comparable with those obtained in 
Refs. [11-13]. 

At this point one comment is to be in order. At first sight it seems that the 
calculation presented in this Letter based on the model with the effective gluon 

7 Notice that fj, > 1.5 GeV is required so that 9 jf dQ 2 T Q 2 T exp(-3Q2) = 1. 
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propagator of the exponential form _D exp(— Q 2 ) can be only used to estimate 
the contribution coming from the relatively low Q T . However, it is not the case. 
The reason is following. As was mentioned before we take for the Tevatron en- 
ergy S'gap/ (G 2 /47r) 2 = 0.6 [10, 16]. If we assume the gap survival factor S 2 ap to be 
0.05 [12, 19] we obtain the non-perturbative coupling G 2 /Atc to be 0.3. From the 
magnitude of the total cross section we conclude that D G 2 = 30 GeV" 2 [2] what 
allows us to extract D = 8 GeV -2 . Now if we compare the numerical values of 
our non-perturbative gluon propagator 8exp(— Q 2 ) with the perturbative one l/Q 2 
we observe that 8exp(-<2 2 ) > l/Q 2 for 0.15 < Q 2 < 3.3 GeV 2 . It means that 
the results presented in Table 1 include contribution coming not only from the non- 
perturbative region, say Q 2 < 1 GeV 2 , but also from the quasi-perturbative region 
1 < Q 2 < 3.3 GeV 2 . Moreover, we have checked that the contribution coming from 
the perturbative region Q 2 > 3.3 GeV 2 with the non-perturbative gluon propagators 
replaced by perturbative ones is negligible. 

At the end, let us notice that the result da/dy(y = 0) = 45 nb for DPE exclusive 
Xc production is consistent with the preliminary CDF Run II upper limit [17], see 
Table 2. Indeed, if we integrate the differential cross section (6) over the rapidity 
range |y| < 0.6 [17] we obtain the value 8 55 nb to be compared with the CDF upper 
limit 80 ± 30(stat) ± 70(syst) nb [17]. 



= 2 TeV 

\y\ < 0.6 


CDF upper limit 


Model 


xl 


80 ± 100 nb 


55 nb 



Table 2: Comparison of the CDF upper limit for DPE exclusive x c meson production 
with the result obtained in the presented model. The satisfactory consistency is 
observed. 



5 Conclusions 

In conclusion, in this paper we investigated the central inclusive and exclusive DPE 
process of Xc> Xb mesons production in pp (pp) collisions. We observed that the 
relation <7 exc i./cc.inci. ~ 0.3 (Xc)> 0-07 between exclusive and central inclusive 
cross sections are hold. The recently proposed exclusive extension of the Bialas- 
Landshoff model was found to be consistent with the preliminary CDF upper limit 
on double diffractive exclusive x° c production cross section. 

8 Since the differential cross section (6) weekly depends on the rapidity of the produced meson 
it is enough to multiply d<r/dy(y = 0) = 45 nb by a factor 1.2. 
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